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Abstract 

The inclusive jet differential cross section has been measured for jet transverse 
energies, Et, from 15 to 440 GeV, in the pseudorapidity region 0.1< |r7| <0.7. The 
results are based on 19.5 pb~^ of data collected by the CDF collaboration at the Fer- 
milab Tevatron collider. The data are compared with QCD predictions for various 
sets of parton distribution functions. The cross section for jets with Et > 200 GeV 
is significantly higher than current predictions based on O(a^) perturbative QCD 
calculations. Various possible explanations for the high-£^T excess are discussed. 
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We present a precise measurement of the inclusive differential cross section for jet 
production in pp collisions at 1.8 TeV. Our measurement is compared to next-to-leading 
order (NLO) perturbative QCD predictions |l|] for jet transverse energies, Et, from 15 to 
440 GeV in the central pseudorapidity region 0.1< It^I <0.7, corresponding at highest Et 
to a distance scale of 0(10^^'^) cm. 

The predictions depend on details of the parton distribution functions (PDFs) and on 
the strong coupling constant as- Our measurement provides precise information about 
both 0, ^. Apart from these theoretical uncertainties, deviations of the predicted cross 
section from experiment could arise from physics beyond the Standard Model. In par- 
ticular, the presence of quark substructure would enhance the cross section at high Et- 
Previous measurements of inclusive jet production were performed with smaller data sets 
by CDF I, § and at lower energy by UA2 and CDF §. 

The measurement described here is based on a data sample of 19.5 pb~^ collected 
in 1992-93 with the CDF detector Q| at the Tevatron collider. The data were collected 
using several triggers with jet Et thresholds of 100, 70, 50 and 20 GeV. The 70, 50 
and 20 GeV triggers were prescaled by 6, 20 and 500, respectively. Cosmic rays and 
accelerator loss backgrounds were removed with cuts on event energy timing and on 
missing transverse energy, as described in reference 0. The remaining backgrounds are 
conservatively estimated to be <0.5% in any Et bin. 

Jets were reconstructed using a cone algorithm[^ with radius R = (Ar/^ + A0^)^/^ = 
0.7. Here r] = — ln[tan(6'/2)], where 6 is the polar angle with respect to the beam line 
and (J) is the azimuthal angle around the beam. The QCD calculation used a similar 
algorithm [|I|. The ambient energy from fragmentation of partons not associated with the 
hard scattering is subtracted. No correction is applied for the energy falling outside the 
cone because this effect is modelled by the NLO QCD calculations. 

The measured jet Et spectrum is corrected for detector and smearing effects caused 
by finite Et resolution with the "unsmearing procedure" described in [0 . A Monte Carlo 
simulation, based on the ISAJET[rU| program and Feynman-Field|rT| jet fragmentation 



tuned to the CDF data, is used to determine detector response functions. A trial true 
(unsmeared) spectrum is smeared with detector effects and compared to the raw data. 
The parameters of the trial spectrum are iterated to obtain the best match between the 
smeared trial spectrum and the raw data. We parameterize the unsmeared inclusive jet 
spectrum with the functional form 

= P„ X (1 - .rf X 10''(-?™>, (1) 

where F{x) = ZlLi -^j ^ pogl^^)]* with E^^^^ in GeV, Po---Pq are fitted parameters and 
xt is defined as 2£'^^"^/ ^/s. The resulting fit of the smeared true spectrum to our data 
yields x^/degree-of- freedom = x^/d-f- = 29.9/34. The best-fit set of parameters for Eq. [|, 
i.e. the "standard curve", are in Table|l]. Corrections to the measured Et and rate for 
each bin of the raw spectrum are derived from the mapping of the standard curve to the 
smeared curve. The corrected cross sections and statistical uncertainties are in Fig.|l| and 
in Table 0. 
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To evaluate systematic uncertainties, the procedure in reference is used. New pa- 
rameter sets for Eq. ^ are derived for ±1 standard deviation shifts in the unsmearing 
function for each source of systematic uncertainty. The parameters for the eight largest 
systematic uncertainties are in Table |I]. They account for the following uncertainties: (a) 
charged hadron response at high Pt'-, (b) the calorimeter response to low- Pt hadrons; (c) 
±1% on the jet energy for the absolute calibration of the calorimeter; (d) jet fragmen- 
tation functions used in the simulation; (e) ±30% on the underlying event energy in a 
jet cone; (f) detector response to electrons and photons and (g) modehng of the detector 
jet energy resolution. An overall normalization uncertainty of ±3.8% was derived from 
the uncertainty in the luminosity measurement (±3.5%) and the efficiency of the accep- 
tance cuts (±1.5%). Additional tests of the unsmearing procedure, including use of the 



HERWIG Monte Carlo program [p!3| to model jet fragmentation, were performed and the 
resulting variations were found to be small. Fig. ^(a-h) shows the percentage change from 
the standard curve as a function of Et for each uncertainty. 

In Fig.|l|the corrected cross section is compared with the NLO QCD prediction |l| using 
MRSDO' PDFs[0,with renormalization/factorization scale = Et/^. These results 
show excellent agreement in shape and in normalization for Ej- < 200 GeV, while the 
cross section falls by six orders of magnitude. Above 200 GeV, the CDF cross section is 
significantly higher than the NLO QCD prediction. These data are consistent with our 
previous measurement 0], which also shows an excess over NLO QCD for the Et > 280 
GeV region. A similar excess is observed when we compare CDF data with HERWIG 
Monte Carlo predictions. 

The distributions of the physical variables in the 1192 events above 200 GeV were 
examined carefully. Data distributions sensitive to the mismeasurement of jet E^, such as 
unbalanced jet Ej- in dijet events, show good agreement with detector simulation. To look 
for time and luminosity dependent variations (instantaneous luminosity increased with 
time), the data were divided into seven time-ordered parts and analyzed independently. 
No significant time dependence was observed. Finally, these events were individually 
scanned and no anomalies were discovered. 

No single experimental source of systematic uncertainty can account for the high-i?T 
excess. For example, in order to reconcile the measured CDF spectrum with NLO QCD 
(MRSDO', /i = Et/2) predictions, we would have to change the jet Et scale by an amount 
ranging from 0.2% at 175 GeV to 5% at 415 GeV, while keeping the change less than 0.1% 
between 50 and 160 GeV. No known feature of the detector, its calibration or the data 
analysis permits such a change. The effects of all possible combinations of the systematic 
uncertainties are included in the comparison described below. 

To analyze the significance of this excess we use four statistical tests: signed 



and unsigned Kolmogorov-Smirnov |14|, Smirnov-Cramer-VonMises |14], and Anderson- 
Darling |TB|, |TB[. For this comparison we choose the MRSDO' PDFs which provide the 
best description of our low Et data. The eight sources of systematic uncertainty are 
treated individually to include the Et dependence of each uncertainty. The effect of finite 
binning and systematic uncertainties are modelled by a Monte Carlo calculation. The 
statistical tests over the full Et range are dominated by the higher precision data at low 
Et] therefore, we test two ranges. Between 40 and 150 GeV, the agreement between data 
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and theory is >80% for all four tests. Above 150 GeV, however, each of the four methods 
yields a probability of 1% that the excess is due to a fluctuation. 

We have considered various sources of uncertainty in the theory. The NLO QCD 
predictions have a weak dependence on the renormalization/factorization scale fi. The 
change in yU scale from 2 Et to Et/4: changes the normalization but maintains the shape for 
Et > 70 GeV [|l^. For the NLO QCD calculations the renormalization and factorization 
scales have been assumed to be equal. Varying these scales independently also has little 



effect on the shape of the theoretical curve [ 18 1. However, soft gluon summation may lead 



to a small increase in the cross section at high E'rlTl, addition, the effect of higher 

order QCD corrections is not known. 

The fractional difference between the MRSDO' [0 NLO QCD predictions and pre- 
dictions using different choices of published PDFs, with = -E't/2, is shown in Fig.[l|. 
The excess of data over theory at high Et remains for CTEQ2M|21|], CTEQ2ML0, 
GRV94|22|, MRSA'H and MRSG||] parton distributions. The variations in QCD pre- 



dictions represent a survey of currently available distributions. They do not represent 
uncertainties associated with data used in deriving the PDFs. Inclusion of our data in 
a global fit with those from other experiments may yield a consistent set of PDFs that 



accommodate the high-ii^r excess within the scope of QCD y, p5 |. 

The presence of quark substructure could appear as an enhancement of the cross 
section at high Et- This effect is conventionally parameterized in terms of a contact 
term of unit strength between left-handed quarks, characterized by the constant Ac with 
units of energy ||2^. While NLO standard model QCD predictions have been available for 
many years, no calculation for compositeness at next-to- leading order [0(a^)] is available. 
Therefore, we have compared our data to a LO QCD calculation including compositeness 
(using MRSDO') and have taken the approach of reference We normalize the predicted 
cross section to the data over the Et range 95-145 GeV, where the effect of the contact 
term with Ac > 1000 GeV is small. The best agreement between this calculation and our 
data above Et > 200 GeV is for Ac = 1600 GeV. This hypothetical contact interaction is 
also expected to lead to dijet production with a more central angular distribution, and this 
analysis is underway. However, until a realistic method for representing the theoretical 
uncertainties from higher order QCD corrections and from the PDFs is found, any claim 
about the presence or absence of new physics is not defensible. 

In summary, we have measured the inclusive jet cross section in the Et range 15-440 
GeV and find it to be in good agreement with NLO QCD predictions for Et < 200 GeV 
using MRSDO' PDFs. Above 200 GeV, the jet cross section is significantly higher than 
the NLO predictions. The data over the full Et range are very precise. They provide 
powerful constraints on QCD, and demand a reevaluation of theoretical predictions and 
uncertainties within and beyond the Standard Model. 
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Table 1: Parameters of the curves corresponding to ± 1-standard deviation changes in 
the systematic uncertainties. 





Po (nb/GcV) 


Pi 


P2 




Pa 


P^ 


P& 


Standard 


3.090x10^ 


-4.128 


1.084 


-0.845 


0.136 


0.00279 


6.733 


High Pt pion(+) 


3.000x10^ 


-4.110 


1.083 


-0.847 


0.135 


0.00213 


6.500 


High Pt pion(-) 


3.118x10^ 


-4.132 


1.084 


-0.844 


0.137 


0.00299 


6.758 


Low Px pion (+) 


3.135x10^ 


-4.163 


1.082 


-0.843 


0.138 


0.00342 


7.209 


Low Pt pion (-) 


3.060x10^ 


-4.096 


1.085 


-0.847 


0.135 


0.00216 


6.272 


1.0% E. scale(+) 


3.174x10^ 


-4.122 


1.083 


-0.846 


0.136 


0.00270 


6.434 


1.0% E. scale(-) 


3.066 xlO^ 


-4.140 


1.084 


-0.844 


0.137 


0.00274 


7.082 


Fragmentation (+) 


3.152x10^ 


-4.161 


1.082 


-0.843 


0.138 


0.00335 


7.214 


Fragmentation (-) 


3.044x10^ 


-4.095 


1.085 


-0.847 


0.135 


0.00220 


6.229 


Underly. Energy(+) 


6.630 xlO« 


-4.314 


1.067 


-0.840 


0.141 


0.00503 


8.045 


Under ly. Energy (-) 


1.730x10^ 


-4.004 


1.099 


-0.846 


0.134 


0.00122 


6.074 


Electron/7(+) 


3.102x10^ 


-4.123 


1.084 


-0.845 


0.136 


0.00271 


6.629 


Electron/7 (-) 


3.106x10^ 


-4.138 


1.083 


-0.844 


0.137 


0.00287 


6.873 


Resolution (+) 


2.422 xlO^ 


-4.082 


1.090 


-0.845 


0.136 


0.00222 


6.645 


RcKolulicni ( ) 


4.262x10''^ 


-4.201 


1.076 


-0.843 


0.138 


0.00306 
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Table 2: The mean true jet Et, cross section and statistical uncertainty. 



<Et> 


Cross Section 


<Et> 




Cross Section 




(GeV) 


(nb/GeV) 




(GeV) 




(nb/GeV) 






14.5 


( 1.14±0.03) 


X 


10^ 


133.8 


( 8 


50±0.12) 


X 


10" 


-•2 


20.3 


( 2.31±0.12) 


X 


10^ 


139.2 


( 6 


62±0.10) 


X 


10- 


-2 


26.9 


( 6.30±0.56) 


X 


10^ 


144.5 


( 5 


00±0.08) 


X 


10" 


-2 


33.3 


( 2.36±0.09) 


X 


10^ 


149.9 


( 3 


92±0.07) 


X 


10" 


-2 


39.5 


( 1.02±0.01) 


X 


10^ 


155.3 


( 3 


13±0.06) 


X 


10" 


-2 


45.5 


( 4.89±0.06) 


X 


10^ 


160.7 


( 2 


46±0.05) 


X 


10- 


-2 


51.3 


( 2.61±0.04) 


X 


lOi 


168.4 


( 1 


75±0.03) 


X 


10- 


-2 


57.0 


( 1.42±0.03) 


X 


lOi 


179.2 


( 1 


10±0.02) 


X 


10- 


-2 


62.7 


( 8.62±0.21) 


X 


10° 


189.0 


( 7 


34±0.20) 


X 


10- 


-3 


68.3 


( 5.43±0.16) 


X 


10° 


200.7 


( 5 


11±0.17) 


X 


10- 


-3 


73.9 


( 3.24±0.13) 


X 


10° 


211.5 


( 3 


41±0.13) 


X 


10- 


-3 


79.4 


( 2.05±0.10) 


X 


10° 


224.6 


( 2 


25±0.09) 


X 


10- 


-3 


85.0 


( 1.44±0.02) 


X 


10° 


240.9 


( 1 


14±0.06) 


X 


10- 


-3 


90.5 


( 1.02±0.02) 


X 


10° 


257.2 


( 6 


67±0.47) 


X 


10- 


-4 


95.9 


( 6.94±0.13) 


X lO'i 


273.5 


( 4 


31±0.38) 


X 


10- 


-4 


101.4 


( 5.18±0.11) 


X 10^1 


292.0 


( 2 


50±0.25) 


X 


10- 


-4 


106.8 


( 3.64±0.05) X 10"^ 


313.7 


( 1 


35±0.19) 


X 


10- 


-4 


112.2 


( 2.64±0.04) X 10-^ 


335.3 


( 6 


37±1.30) 


X 


10- 


-5 


117.6 


( 2.00±0.04) X 10^^ 


364.0 


( 3 


03±0.66) 


X 


10- 


-5 


123.0 


( 1.48±0.03) X 10-1 


414.9 


( 9 


05±2.86) 


X 


10- 


-6 


128.4 


( 1.10±0.03) X 10-1 
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Figure 1: The percent difference between tlie CDF inciusive jet cross section (points) 
and a next-to-leading order (NLO) QCD prediction using MRSDO' PDFs. Tfie CDF 
data (points) are compared directly to the NLO QCD prediction (hne) in the inset. The 

normalization shown is absolute. The hatched region at the bottom shows the quadratic 
sum of correlated systematic uncertainties. NLO QCD predictions using different PDFs 
are also compared with the one using MRSDO'. 
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(b) Low P,j Hadron response 



(d) Fragmentation 
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Figure 2: The percentage change in the inclusive jet cross section when various sources of 
systematic uncertainty are changed by ±l-standard deviation from their nominal values. 
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